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Abstract 

A  technique  is  discussed  by  which  asymptotic  solutions  of  re- 
flection problems  in  paraboloidal  and  parabolic  cylinder  coordinates 
can  be  found  for  large  wave  numbers.  This  technique  applies  to  re- 
flections from  the  interior  of  parabolic  and  paraboloidal  reflectors. 

Detailed  discussions  are  then  given  for  the  cases  of  reflection 
of  incoming  plane  waves,  waves  emitted  by  sources  at  the  focus,  and 
waves  emitted  by  sources  on  the  axis  of  the  reflector,  for  both  co- 
ordinate systems.  For  the  last  problem  the  nature  of  the  caustic 
surface  and  the  behavior  of  the  reflected  wave  on  the  caustic,  off 
the  caustic,  and  at  the  cusp  of  the  caustic  are  discussed. 
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1.     Introduction 

The  purpose  of  this  report  is  to  discuss  the  asymptotic  solutions  of  some 
diffraction  problems  in  parabolic  and  paraboloidal  coordinate  systems,  in  the 
sense  that  the  wave  nunber  k  will  be  large, 

Buchholz*-  ^  has  derived  Green's  function  representations  in  these 
coordinate  S3rsteins.  Using  his  formulas  it  is  possible  to  derive  the  exact 
solution  to  the  problem  of  the  reflection  of  a  spherical  wave  on  the  inside  of 
a  paraboloidal  mirror.  This  leads  to  a  contour  integral,  which  does  not  lend 
itself  readily  to  asymptotic  expansions.  In  the  neighborhood  of  the  cusp  of  the 
caustic  surface  the  behavior  of  the  reflected  wave  is  such  that  the  solution 
depends  very  stixjngly  on  how  the  cusp  is  approached.  The  exact  solution  indicates 
only  one  of  several  possible  patterns  of  behavior. 

Various  techniques  have  been  developed  to  find  the  asymptotic  solutions 

r2i 

of  reflection  problems.  Keller,  Lewis  and  Seckler"-  •'  have  solved  numerous 
problems,  but  not  where  caustic  surfaces  were  involved.  An  approach  to  the 
asymptotic  solution  of  diffraction  problems  in  two  dimensions  where  caustic 
surfaces  are  involved  was  discussed  by  Kay  and  Keller*-  •'. 

The  technique  developed  in  the  present  report  allows  one  to  construct 
asymptotic  solutions  of  the  wave  equation 

^u  +  k  u  -  0 

where  the  bo\indary  values,  u,  are  given  on  the  inside  surface  of  a  paraboloid  of 
revolution  or  a  parabolic  cylinder.  The  solution  can  be  expressed  by  Green's 
theorem  by  the  integral 

^-  dG  , 
u^ds. 
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Asyiaptotic  expansions  of  -^r:  are  used  in  the  integrand  to  construct  the  solution, 

Oli 

Th£  resulting  integrals  can  be  evaluated  readily  when  points  of  stationary  phase 
can  be  found. 

A  similar  technique  is  developed  for  probleins  where  ttj  is  prescribed  on 
the  boundary. 

It  should  be  observed  that  these  asymptotic  forms  show  that  Huygens'  principle 
has  scrnie  solidity  in  these  cases. 

The  following  problems  are  treated: 

a)  reflection  of  incoming  plane  waves  from  the  interior  of  »  paraboloid 
of  revolution  and  a  parabolic  cylinder; 

b)  reflection  of  spherical  waves  due  to  a  point  source  at  the  focus  of  a 
paraboloid  of  revolution  ,  and  reflection  of  cylindrical  waves  due  to  a  line  source 
at  the  focal  line  of  a  parabolic  cylinder; 

c)  problems  similar  to  b)  except  that  the  sources  are  on  the  axis  of  the 
reflector  away  from  the  focus  or  focal  line. 

In  the  last  case  different  asjinptotic  expressions  are  obtained  for  points 
on  and  off  the  caustic  and  at  the  cusp  of  the  caustic. 
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2,   The  caustic  surface 

If  one  considers  the  case  of  a  point  source  placed  on  the  axis  of  the  re- 
flector one  finds  that  the  envelope  of  the  reflected  rays  plays  an  important  role. 
That  envelope  is  called  the  caustic  surface. 

We  will  now  derive  the  equation  of  the  caustic  surface  and  derive  some  of 
its  properties.  Consider  a  (p ,z)   coordinate  system  with  a  parabolic  reflector 
given  by  the  equation 

2 

where  y>     denotes  the  focal  length  and  the  focus  is  placed  at  the  origin.  A  point 
source  is  now  placed  at  the  point  (0,a)  and  it  can  be  shown  by  geometrical  optics 
that  the  ray  which  intersects  the  reflector  at  (p,z)  will  be  reflected  into  the 
surface  and  will  intersect  the  axis  at  (0,q),  where 

q  .  i ^  ''   . 

The  slope  of  the  reflected  rsy   is  given  by 

q  -  * 

P  -  ^  _    . 

-P 
Therefore  the  equation  of  the  reflected  ray  is  given  by 

z  ■  q  +  pp  . 

This  line  must  be  tangent  to  the  caustic  surface  at  some  point  (p  ,z  ),  so  that 

t^     .     q  >  p^p 

where  z  -  f(p)  represents  the  caustic  surface.  Since  one  can  asscciate  one  point 

(PjZ)  with  every  (p  ,z  )  we  can  write 

c  c 


-  u  - 


x>,/      \  dq  dz       dp  dz 

f  »(.o  )  -  _a._+  P„—  —  +P  -  P, 


C        ■"  c   ^ 

dz  dp       dz  dp 


from  which  we  see  that 


C        p'   * 


p'q  -  q'p 
*c     V' ' 

and  since  q  and  p  are  functions  of  z"  thase  represent  the  parametric  equations  of 
the  caustic  surface. 

Those  points  where  z  5j  -  y?    correspond  to  rays  which  are  reflected  from  points 
near  the  vertex  of  the  parabola.    One  can  now  show  that  for 

z  »   -  y|  , 

P    ?i    0 

the  caustic  surface  is  given  by 


p-5 

p 

v^(a  +v^) 

2       -2 

a 

?'■ 

U>1  (z  +  v^  ) 

Now  p  can  be  eliminated  and  we  obtain 

which  shows  that  at  the  point  (0,  -* —  )  the  caustic  has  a  cusp.  Figure  1  shows 
the  general  nature  of  the  caustic  surface. 
We  now  consider  the  function 
R  *  r, 
where 


-  Aa  - 


(0,-^) 


Flgyire  1 
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V  (z  -a) 


5 — =r 

♦  P 


r  -  y'  (z  -z)  +  (p  -  py 


and  seek  those  values  of  z   for  which  it  is  stationary.  We  will  assume  that  p  is 
negative,  (p,z)  lies  above  the  caustic  ,and  p  is  positive.     Then  R  +  r  must 
have  three  stationary  points.  Two  of  these  can  be  located  by  geometrical  optics, 
because  they  correspond  to  paths  of  least  time.  If  it  is  possible  to  draw  tangents 
from  (p,z)  to  the  caustic,  then  the  intersections  of  the  tangents  with  the  para- 
bolic surface  correspond  to  paths  of  least  time.  But  since  ^   "^  is  negative  at 
z»-")0,p«Oit  follows  that  R  +  r  Btust  have  a  local  maximum  for  some  point  where 
p  is  less  than  0  but  greater  than  the  point  corresponding  to  the  path  of  least  time. 
If  p  ■  0  then  the  local  maximum  corresponds  to  the  point  (O,  ->?  ).  In  this  case  the 
stationary  points  are  symmetrically  located  on  the  parabola. 

As  the  point  (p,z)  approaches  the  caustic,  the  two  stationary  points  for  p  <  0 

approach  each  other.  Then  in  the  limit  they  merge,  so  that  not  only   d(R+r)/dp 

2      -2 
but  also  d  (R+r)/dp  vanishes.  The  stationary  point  is  a  point  of  inflection, 

so  that  d-^(R+r)/dp-^  does  not  vanish. 

For  points  below  the  caustic  only  one  stationary  point  is  left  for  which 

p  >  0.  For  points  on  the  axis  and  below  the  caustic  that  point  is  (O,  -  ^  )♦ 

3,   The  Green's  function  in  parabolic  coordinates 

The  parabolic  coordinate  system  can  be  defined  by  the  transformation 
X  •  2V  gr^   cos  9 
y  »  2  VTt"  sin  J8 
z    '    4  -  ^  . 
The  surface  y|  ■  y|   defines  a  paraboloid  whose  axis  coincides  with  the  z- 
axis,  and  which  includes  the  entire  positive  part  of  the  z-axis.  Similarly  C  "  4j^ 

is  also  parabolic,  but  opens  in  the  opposite  dirsction.  When  expressed  in  parabolic 

2 
coordinates,  the  wave  equation  ^u  +  k  u  "  0   beconies 
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.2 
+     k     u 


0   • 


Separation  of  variables  then  leads  to  the  system  of  ordinary  differential 
equations 


dfAK) 


3C 


K    -V-    *    (""^'k   '  ''*^^)  'l^^^    ■    ° 


f,    >(^^     (.^.^^*2i.^)f,(v,).    0 


X-yy      *        H'     f     (I?)        -       0, 

d9  ^ 


where  ]C   and  (ji  are  arbitrary  complex  parameters.  In  the  notation  of  Buchhola'»  J, 
the  two  linearly  independent  solutions  of  the  first  of  these  are: 


t^{K)   -  ■'J^  (-2ik4)  - 


r  (1  *  ti) 


where  the  Kumraer  function  is  defined  as  usual  by 

00       (a)^    ^n 
^F^(ajb,z)     -     £     ^    ^ 
n"0  n 


(a)       ■     a(a+l)    ...   (a+n-1) 
n 


(a). 


and 


f^CC) 


w^  (-2ik4) 


-  .'^^(-2ik4) 


sinnp. 


.J(-2ik^) 


_r(^./)  "  n^-iL), 


When  n  is  an  integer  the  function  \it{z)   must  be  derived  from  the  above 
definition  by  a  limiting  process.  Similarly  the  two  linearly  independent  solu- 
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tlons  of  the  second  equation  can  be  stated  as 


fgCv^)  -  ia^2ik>|) 


(2ikv^)'^/2a-^^^  ^Fi(  ijii  -  X  Jl*w2ik  >^) 


r  (1  ♦  t^) 


and 


Vn)-<(2ik>l)-g^ 


f  m-^^'(2ikvi)  H^(2ikn) 

In^'ic)  '   rciji - /^  )  J' 


When  M.  is  an  integer  iiC(z)  is  regular  and  single-valued  over  the  whole  space, 
but  w^z)  is  in  general  neither  single-valued  rwr  regular,     Buchholz  ^  -I  has 
shovm  that  the  Green's  function  for  the  first  boundary- value  problem  has  the 
following  integral  representation  in  parabolic  coordinates! 


P-<»  i,r-ioo 

p 


l+Px 


(8.  i^)r(-8-  ^) 


P(.2iky^^)     wP(-2iky^^) 
|P(-2ikv,^)     w^(-2iky|^) 


m 


P(-2ikr,    )         inP(-2ikC)       w^  (-2ik£   )    1  K<1 

s  /     -s  "^  °      I  ds  , 


mP(-2ik>)    )         mP  (-2ik^^)     w^  (-2ik4) 
s  I  a        I      — s  w         — B 


=0 

'o 


^ai     ^0=^  ^    ^°'  0<"^p<T^- 


By  use  of  the  Wronskian 


&] 


W^  wj(z),  bP(z) 


rp^-s) 


it  is  possible  to  widte 
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30^ 


^K  ''^. 


.2ik     ^        iP^^'-^^o^   1 

5m7  ^  „®  ?;n: 


-(f'+iOO  .    ^ 


la    p«-oo 


V-ioo 
P 


mP(.2ikY|     ) 


mP^(-2ik^)  wP„(-2ikt   ) 


-s 


-8 


e      <       4, 


mP  (-2ik£   )wP  (-2ik4) 


\    ds 


^    >     ^o   • 


The  integrals  under  the  summation  can  be  evaluated  by  closing  the  contour 
on  the  left  and  suimning  over  the  residues  of   P(s  +  -5^).     Heference  to  the 
asymptotic  fonmilaB '-  -»   for  the  parabolic  functions  with  respect  to  the  parameter 
s     shows  that  this  is  permissible  provided  that 

This  leads  to 


9G^ 


%-\ 


sir  t  y'''""'  t  ^^^"^ ■!.^<-^'-L W^^^>" 


la     p«-co 


n«o 


^ 


n+ 


n+ — - 


n+— ij*- 


It  now  follows  from  the  definition  of  the  function  in^(z)  that 

7- 


B-   ,_(-2ik£  ) 


n+ 


1+p^ 


-P/2 


-"^     W^ii^^a) 


-^oV   i^^Vna^  ^Fi(-n;l.p|-2ik4^) 
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The   KuDuner  functions  appearing  in  the  fraction  are  polynoraiais  of  degree  n. 
If  it  is  now  assuBjed  that  k^     and  k>7      are  large,   one  can  expand  this  ratio  in  an 
asymptotic  series  in  k^     and  kvj    .     We  assume     for  the  moment,  however,  that  al- 
though k  becomes  large,  ^  and    Y)     are   such  that  k^  and  k  v?     can  still  be  consider- 
ed relatively  small.     Then 

^F^(-n}l+pj-2ik4^)      ^/     C^\n 
^F^(-nil+Pj2iki^^)      "  \^        a 

Now  use  of  the  generating  function L  J 

n*o  n+  -—■       n+  -it-  v'      \±+t) 

allows  one  to  write  asymptotically 


9G 


exp  I  ik 


iv'?.'*(«-Vi7ir]''. 


V«o  1) .  ('*  V^T^ 


o  '  a 


PV   ^o  *  ^ 


a  ^o  -iy   \  ^0    'a 


*a 
This  expression  can  be  simplified  by  use  of  the  generating  function  of  the 

Bessel  function,  i.e., 

00 

p—oo 
then  one  obtains 


YZ        J  (z)tP  -  exp|z(t  -  ^)  /  2]   ; 


30^ 


^'^'>.^o'^. 


,_-^e.(..  [,..V(C-V  ^-  ^-^co.  <.^„,] 
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The  distance  between  two  points  in  parabolic  coordinates  is  given  by  the  ex- 
pression 

r  •/(^a*^o^^  *  (n+«>^  +  2(C->|)(>|3-4o)  -  8V^y|C^^;  cos  (g(-(?^) 

which^ under  the  assumption  that  the  quantities  "^  _»  4  are  large  relative  to 

v|  ,   ^^can  be  expanded  in  a  binomial  series.  Use  ol  w^at  series  shows  that  one 

can  now  write 


■^o""i 


•2ik  e__ 
m    r 


ikr 


o  >a 


and  inserting  this  result  in  the  expression 


u  ■ 


3GjL 


ds 


leads  to 


2  v^ 


-  90  "1 


^o-^< 


'o   o 


In  a  similar  manner  it  is  possible  to  start  with  the  integral  representation 
of  the  Green's  function  for  the  second  boundary -value  problem  to  derive  an  asymp- 
totic representation.  Then  it  follows  from 

-<Jl  +  ioo 

.    -2ik  ^    ^p(^-<^o>  1    r 


e 


P--00 


5^  j     r(8*^)r(-s+^)»>e(-2ijc4)wj(-2ikc^) 


-C-ioo 
P 


"^s^-^^^o^    ■!«(-2ikv|^) 


m^  (.2ikr|) 


ds 


»rj-2ik  y^j 


that  asymptotically 


-  il  - 


and 


ikr 


u 


1   e__ 

?n   r 


n 


.00  2n   -I 

i  if 

^0   o 


ikr 
—-dCdgf 


U»   The  Green 'e  function  in  parabolic  cylinder  coordinates 

The  parabolic  cylinder  coordinate  system  can  be  defined  by  the  transformation 


X     • 


4->7 


y    =    ±    2  /lY 


p    -    yx  +y 


CM 


X     ■     p  cos  ^  y     ■     P  sin  0  • 

As  in  the  previous  case  the  surfaces  ^  ■  constant  and  ^  «  constant  are  parabolas. 
The  ^ave  equation 

2 

^u  +  k  u  »  0 

when  expressed  in  parabolic  cylinder  coordinates  becomes 


iK 


^K 


'^^^^  )  7^  *  mTtF 


+  k  u  ■  0  • 


Separation  of  variables  leads  to  the  two  ordinary  differential  equations 
f'l  (/!")     +      [(2k  v^  )2  +  Uik(Y  *  ^)]  fi(  Vr  )     -     0 

fg   (V^)     +      [^2k  v^  )2  -  Uik(Y  +  5-)]  f2(  V^  )     -     0, 
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where  y  is  an  arbitrary  complex  parameter.  Two  linearly  independent  solutions 
to  the  first  of  these  can  be  defined  by 


E 


^°^  (2V::ik^)  -  V57  (-2ikO^/'*  m:^/\  (-2ik4), 


i*  IT 

Another  solution  which  proves  convenient  can  be  expressed  in  tenns  of  these  two 
as 

D^(2  /IHE^  )  -  2^/2  ^^^     1  JC — ^ J^ 


-  2Y/2(.2ik4)^/^  wy^  A-2±kK), 

where  the  functions  »y(z)  and  Wv(z)  were  defined  in  the  previous  section.     Similarly 
the  solutions  of  the  second  equation  can  be  written  as 


E^°\2VlEyj)  -y2ir  (2ikV|)l/'^m;^/\(2ikvj), 


-  2^/^(2ikV|  )^/'*  w^/^  l(2ikv|  )    . 


-  13  - 

Buchholz  L  J  has  shown  that  it  is  possible  to  represent  the  Hankel  functions  of 
the  first  as  well  as  the  second  kind  in  terms  of  these  functions.     The  former  can 
be  expressed  by 


-ni 


H^^^(kr) 


/2e 


1 


-6+ioo 

I 


r(-8+^)r(s*^) 


E^°^    ^(2/1114   )   D  ^(2/^11^); 

-2s-j 

(0) 


.2S-J 


-^-ioo 

|<r|  <  i/u 


E^^'      (2/-ik£    )D 
■2s-i. 


-.»4 


(2/:ik4  ) 


E^°\(2/r3[F7-)     D^    ^(2/^Ik7i^  ) 

2s- ^ 


28-5 


ds 


.(0) 


E^'''',(2y-ikv|  )  D     .(2/:iir?j ) 


2s- 1 


2s-: 


^  -6+ioo 

-t!i  / 

^2e  1 

"n  ?ii 


r(-s*  J)r(s*  I) 


-6'-ioo 
1^1  <  3/U 


E^-^^  ,(2/:Tkr)    D      i(2/::ikr-) 

-2s-^  -26-| 


-2s.| 


.2s- 1 


.(1) 


E^-"^    ,(2/^Ikii)       r         .(2K-ikv|^  ) 
+28-i-  +28-2- 

E^-^^  ^(2i/::ik^)   D      ^(2/^   ) 


\ 


♦2s- i 


+28-i 


ds 


^  <  ^o>       ^  <    ^ 


^>^o'       V    ^, 


With  this  integral  representation  it  is  now  possible  to  construct  the  Green's 
functions  for  the  first  as  well  as  ths   second  boundary-value  problems.     We  obtain 


-  lU  - 


in  the  case  of  the  first  boundary- value  problem 
-fi+ioo 


^2e""^       1 


T^I 


(-S+  j.)n(s+  \)  E^°^  i(2y:ik4 )  D    ^i^v^m; ) 


-<f-ioo 
|<r|<  lA 

25-2- 


,(0) 


2s- 


-2e- 


D    ^(2  v^iTvf;; ) 


28- 


E^^\(2  V=ir?|7  )  D       ^(2  V^Iik^  ) 


2s>; 


-2s- 1 


E^°\(2  yiiFi ) 

2s4__ 

E^°^,(2VliF^) 
2s-^  ^ 


ds 


-d+ioo 


*^^^    5^     (  r(-s*l)riB*l)E^^\{2V^m)T)         ^(2yiikC^) 


-2s- i 


-6'-ico 

Kl<  3A 


-2s-i 


E^^^  ^(2  y^urr;; ) 


+2s-^ 


.(1) 


D         ^(2  /-ikTi7  ) 


+28-^ 


E^-^'    ^(2  \/-irv|^  ) 
+2s-^ 


D         ^(2  /^Im;  ) 


+2s- 


E^^^    ,(2V^iF|  ) 
+2s-i- 

TO 


ds, 


E^-^^  j(2/:iki[;;) 


+2s- 


a<4^,  i<  n^). 


Use  of  the  Wronskian 


[1] 


w<!  yz),  E^°^(2) 


y^ 


wJ   Djz), 


.    .,e('^(z) 
Y  Y 


■  yiT 


1*    Y/2 
2 

1*y/2 
2 

r(^) 
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allows  one  to  calculate 


-6+100 


Jv^ 


V^    ^    j         2r(-X°^    i(2/:ii^)D         ,(2/11^) 


-(T-ioo 
Iff  1<L/U 


E^(2yrnnr-) 

2s-=. 


d8 


-6+ loo 

375"  /■ 

'■a  ^  -2E-X-  -4:8- «• 

-<r-ioo  ^ 

l<y|<  3/u 


28-5- 


ds 


(€  <  ^o'. 


Whsoaver  we  have 


/TI  >  /«  *  v«. 


1 


we  can  close  the  contour  on  the  left  and  evduate  the  above  integrals  by  susoning 
over  the  residues  of  the  gamma  functions. 


30, 


^< 


oJV^ 


/^ 


i;;?7?  f  00  2-^-ir(.)N  (,) 


,(0) 


o     <  a 


^a     j  fe) 


Nl 


4n  (2  V^^  ^  V2  v^=^  )  ;To7 


El2i.i(2  /^ICT  ) 


^:2N.i(2  vuft;) 


00     „"  'JT/    ^N 


,(1) 


♦  i;^      Ml     ■   niii^^  V:^  ">2N»i'^  ^'^^ 
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B7  ^»se  of  the  functions  in^(z)  and  "^^i^)  it  is  now  possible  to  rewrite  this  ex- 
pression in  a  form  similar  to  the  one  arrived  at  in  the  case  of  the  paraboloid 
of  revolutloni 


vn 


1/U 


n) 


in-^f(.2ikf  ) 
\i'(2iMa^ 


00  r(N+ x)  ,/9        ,.9 


in^/2   (.2ik£   ) 


3(2ikv|    ) 


If  it  is  now  assumed  that  k?  and  kV|   are  large  relative  to  k^  and  k  vi   one  can 
take  advantage  of  the  asymptotic  relations 


±172" 


m-  ■    i(2ik^  ) 


\ 


and  then  by  an  application  of  the  generating  function  of  the  confluent  hypergeo- 
metric  functions  ^^^   one  obtains 


90^ 


Jv^ 


o  ^a 


7a*^o 


exp 


^.<r 


i^(^a*^o*(^->()>^] 


'*ni/U,  f^^j/^Ml 


^V  ^a*^c 


'    '^1    y)    *  K — 
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"Hie  distance  r  between  two  points  can  be  expanded  into  a  rapidly  converging  bi- 
nomial series  if  f     and   r>     are  large  compared  to  g  and  >7     : 


r  -V(^a*^o)^  +  (^   H)^  *  2(4-v|  )i^^<^}  '  8  yTv^ 


>^      .  .     *  U  M   )  ir-TT^    -  _mVLa 


'a      ^o 


Hence  one  can  replace  the  preceding  expression  by 


3^ 


^o-"^, 


-/ 


T5IF  ikr 
— — —  e 
nr 


and  since  this  is  an  asymptotic  expression  it  could  be  written  as 


30^ 


^o-^. 


2ik  H^-"-^  (kr) 


One  can  now  write  by  means  of  Green's  theorem 


-  ^1 


^"la-psr 


ds 


1 

HI 


_  ^^1 

"5T 


'o    ^a 


/^d^-2/ 


^:hfT 


.00 


a 


u    • 


ikr 


/-oo 


V^ 


dyT 


The  above  integrals  must  be  taken  over  the  whole  reflector,  which  falls  into  two 
parts  defined  by 

s/y^  •  /rp  ,   V^  >  0  and  i/C  <  0, 


A  process  similar  to  the  one  Just  carried  out  leads  to  an  asymptotic  repre- 
sentation of  the  Green's  function  for  boundary-value  problems  of  the  second  kind. 
Then 
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■A 


'2e 


^ir 


-6+ioo 


1 
21a: 


r(-s*  ^)  r(s* })  E^°^  ^(2  v=:te4  )d    ^(2  y=iE^ ) 


-28- 


2s-s-  2a-x 


-28- i 


E(°\(2VlIinr)      D       ,(2ylTF^) 
2s-i-  '°  28-i.  '° 


E^°\(2yCiiFf ) 

2s-j 


ds 


-d+icx) 


;  -28-5-  -28-5- 


-5'-ioo 

Kl  <  3A 


£^^^1^(2  VhCF^) 


28-1 


28-2- 


D       ,(2  >/'ikri      ) 
28-^ 


D     3^(2  y-ikv|^  ) 


2s- 


E^^^'(2v:drFT 


ds. 


(«<  V      "^^   ^^o^' 


from  which  one  can  show  by  a  similar  procedure  that  asymptotically 

ik(r  -  f) 


0,  -  2 


yiSor" 


-V?  4^^kr)  . 


Then  it  follows  that 


i  (  _  au  ..     1 


/   ik(r-  ^)      /-— 

yS5^        ^^   f    ^ 


dC 


-i 
ie  ^  ,  ^^.  3^ 


a 


/   ikr  ^- 


'-00 
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5 •       Applications  of  the  asymptotic  forroulaB  to  diffraction  problems 

We  will  now  apply  the  previously  derived  formulas  to  nvunerous     cases. 

a).     First  we  shall  treat  the  diffraction  of  a  plane  wave  by  a  paraboloidal 

reflector.      Then,  taking 

-ik(C  -yj  ^) 


u 
we  have 


ikv,        ,■="   .2.      Ik  (?-«•>?  J 


"  ■  *  nr-^  " — R "«*. 

0         0 

where 


so  that  u  is  evaluated  at  the  point 


y     •     0 

X     =     Co  -  >] 

and  the  general  point  on  the  paraboloid  is  given  by 


y^  "  2  y?  ^  a  sin  9^ 


^a  ■  ^-  ^a- 

Two  points  of  stationary  phase  can  be  found.  Both  lie  in  the  plane  y 
and  are  given  by  the  intersections  of  the  parabola  and  the  line 

!i     X 

which  leads  to  the  quadratic  equation  in  4 
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Danote  the  root  corresponding  to  (Jf  ■  ji  by  4,   and  the  other  by  4_,     Then  near  these 
roots  one  can  use  the  expansions 


p.p.    ^iiVySlir    i^    ^2  ^v^S[vrT 


near  C-i  ^^^cl 


<«-«2'*  /vTTr    .2  2/iivn" 


near  Cp*     It  is  now  possible  to  write  the  solution  by  an  application  of  the  nethod 
of  stationary  phase  •     We  have 


u 


-^       ik(V4iMa)        nTz      ^^\<2'\^^ 


^o^  J     M 


/V2 


The  above  solution  is  valid  only  if  the  coordinate  x  does  not  vanish.  If  it 
does  there  is  no  true  point  of  stationary  phase.  However,  one  can  still  evaluate 
the  integral  asymptotically  in  this  case.     Then 


is  independent  of  (jl^  so  that 

u    -     -2ik  v|  ^  e  *  5-^ d?     . 

o 

If  we  let  X  -  k^  -onder  the  inte^ral^  use  binomial  series^  and  drop  higher-order 
terms  we  obtain  directly 

ik(2v^3*?^) 


e 
u     ■     — 


«o 
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This  analysis  shows  that  the  diffracted  wave  has  the  character  of  a  point 
source  near  the  focus. 

An  analogous  procedure  can  be  followed  with  the  parabolic  cylinder,  and 
one  obtains 

for  off -axis  points.  On  the  axis  for  ">^  ■  0  one  obtains 


u  -  2i  y  - —   e 

/  ^o 

b).  We  will  now  consider  the  problem  of  a  source  at  the  focus  of  the  re- 
flector which  results  in  an  outgoing  plane  wave.  This  problem  can  be  considered 
as  a  reciprocal  problem  to  the  problem  considered  above,  where  an  incoming  plane 
wave  results  in  a  singularity  at  the  focus.  One  proceeds,  as  before,  with  the 


f  ona'il  a 


where 


ikr 

e 


y(^o-  ^  '  ^   Ma^  *  ^^o\    -^  '^^^a  -  ^  ^^^oX^  ^  cos  0  , 
r  -  4  M 


• 


The  point  at  which  the  expression  R  +  r  is  stationary  can  be  determined  by  geo- 
metrical optics  and  is  given  by 

?  -  C   ^ 
°   ^a 

Qf  »  0  . 
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Expanding  R  ♦  r  about  that  point  leads  to  the  following  expression  for  the 
solution  : 


''(na-T>''la*«o'("la*«on 


a  'o  V  '     '-00 


a     I  ' ^  \  a  '0 

-00 


u     -         P^  exp(ik|   2v|^.?^.v^]), 


^      +  C  >1 
^a       ^o  \ 

which  is  similar  in  character  to  an  outgoing  plane  wave,  but  is  decreasing  in 
amplitude . 

An  analogous  operation  with  the  integral 


u 


leads  to  the  function 


r 

ik(R+r) 

6 

n 

ioo 

ySF 

dVC 


-ij  ik(2V|   +Co-V]) 

2  V2^  e    e 
u  ■  — t 


nVk(v^^+  4-) 


^ 


which  is  the  corresponding  solution  for  a  line  source  placed  at  the  focus  of  a 

parabolic  cylinder. 

c).  In  this  section  we  shall  consider  the  case  where  a  point  source  is 

placed  on  the  axis  of  the  reflector,  not  at  the  focus,  but  at  a  distance  from 

it.  As  before,  we  take  as  our  starting  point  the  formula 

Ik  V7     ("^    .2n    ^^^ 
u  -  -  -^        u  V-  ^«  ^> 

4    i 
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where 


—  e 

u     -     -_- 


R     -   /(C  -  >7  ^  -  a)^  *  hi^ 


a     ) 


r    -    /(C-  ^  a"^o+  "^f  *  hK^l^*  UKy]     -  8  V^^^  >|  v|  ^    cos  0 
The  point  given  by 

X    -     2  y?;^  ,        y     -    0,  z     .     4^  -  v^ 

is  the  point  under  examination} 


X  -  2  74  v^  3  cos  gf,  y  -  2  VT^^   sin  (3,    z  -  C  -  v^^^ 

is  a  general  point  on  the  surface  of  the  paraboloid;  and  the  source  is  located  at 

X  ■  0,     y  ■  0,     z  »  a. 

Unless  the  point  under  examination  lies  on  the  axis  of  the  paraboloid  or 
the  caustic  surface,  the  function 

u  -  0(1). 

This  follows  from  simple  stationary  phase  considerations  and  the  fact  that  R  +  r 
Is  such  that 

jVnL  /o  is^  f  0 


on  the  surface,  but  at  some  points 
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simultaneously. 

If  "^    "0,  the  point  under  examination  lies  on  the  axis.     In  this  case 
the  (f  integration  is  trivial  and  our  formula  ireduces  to 

f°°     ik(R+r) 

Jo 
The  expression  R  +  r  becomes  stationary  at 

2 

But  I"  >  0  if  and  only  if  ^  >  v?  /a,  so  that  we  have  a  point  of  stationary  phase 

only  if  the  point  under  inspection  lies  above  the  cusp  of  the  caustic.  Then  the 
integral  can  be  evaluated  by  stationary  phase  techniques  and  one  obtains 


o  /-rt —  ik(R+r) 
u  ■  -  


vt^^ + 4/yF 


Here  IT  and  r  represent  the  values  of  R  and  r  at  the  poliit  of  stationary  phase. 
It  may  be  observed  that  u  is  on  the  order  of  yf  and  remains  finite  even  at  the 
cusp  of  the  caustic.  For  large  values  of  a  the  formula  reduces  to 


-^/? 


ik(a+2V|^) 

e 


at  the  cusp  of  the  caustic. 

2 
For  4  <  >7  /a  a  different  technique  must  be  employed  to  evaluate  the 
0   I  a 

integral,  since  no  point  of  stationary  phase  exists.  We  replace  K  ^V  x/k  under 
the  integral  and  simplify  the  integrand  by  dropping  all  terms  of  order  1/k  and 
higher.  Then  the  integral  can  be  evaluated  explicitly  and  vre  obtain 

^(2>^a*  ^*     ^o^ 

-^a« 
u  ■   — 


la  -  «^o 
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Here   u  is  of  order  i  with  respect  tc  k  and  hss  a  simple  pole  at  the  cusp  of  the 
caustic,  which  is  completely  different  frora  the  situation  when  the  cusp  is 
approached  from  the  opposite  side. 

rtill  a  different  kind  of  behavior  is  observed  on  the  caustic  and  at  the 
cusp  when  the  cusp  is  approached  along  the  caustic. 

If  we  now  choose  £  and  ^  so  that  they  lie  on  the  caustic  surface  we  have 
o      1  " 


where 


u  -  ^   I   j   2— gp dCd0^ 


'o   0 


R  -  >/a-V|^-a)2.U4v^^  , 


>^(^-^a  -  ^o  •'^  ^^  *   ^^na  *   ^^o"!  -  ^  ^^^^TtTT    °°^  ^  ' 


9(R>r)  .  ^/^^oM" 
~9?    "       r 


sin  Q(  -  0 


if  P  -  0,  n. 


For  0  ■  0  and  n  we  can  also  find  points  of  stationary  phase  such  that 

9(R^^)  -  0 

but  for  0  ■  Ti  it  is  also  true  that 

aVrO  .  0 

34 

because  4  and  \    are  on  the  caustic.  Hence  the  only  point  to  be  considered  is 
the  one  for  which  0  »  n.  Then  the  integral  can  be  simplified  to 


ikv^    ,^*^^ik(K.r)     r 

'b  -oo 


-2ik 


y^VJVa^2 


dgr 


^  ''"la  V^ 


°°     ik(K+r) 
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Near  the  point  of  stationary  phase  ^  one  can  expand  R  +  r  in  a  Taylor  series 

R  +  r     -     (R  ♦  r)     ♦     ^^5^    A^  ♦   ... 


where 


P     .     3^(R+r) 


34^      -■  €  -  4 


so  that  we  obtain 


u     -     ^^^  ^         j       exp|ikl-^Jd4 


/-co 


WU^    6-^/3  p(l)^ik(R.r)^l/6 


e 

a 


In  the  immediate  neighborhood  of  the  cusp  of  the  caustic  we  can  use  the 
followiiig  relations  to  simplify  the  expression  for  ul 


ff  ;^  a  +  V| 


a 


r   ':^  Y^ 

a^     <^- 

^0    ^ 

^a/^ 

161"^  a 

\     OC 

(a.n^)'  ^a 

6a 

/..   ..     ^Z    ^ 

Then  one  obtains  for  u  near  the  cusp  and  on  the  caustic 
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u  ■ 


Hence  it  may  be  observed  that  u  increases  with  k  like  k  '  and  that  u  has  a 
singularity  at  the  cusp,  but  yet  the  behavior  is  different  from  that  trtiich  i6 
observed  when  the  cusp  is  approached  from  abo-vo  or  below  along  the  axis. 

To  analyze  the  case  of  a  line  source  on  the  axis  of  a  parabolic  cylinder 
we  proceed  in  a  similar  manner  by  starting  with 


u  -  2  e' 


where 


so  that 


u-  -H^^^(kR) 


> 

r"          ikr 

u   ? —  d  vT 

-co 

^  -r-. 

2-      i(kR-J) 
Co      e 

°°   .ik(R+r) 


*    '        ysr 


u   «   -  /  >? 

''-00 


When  the  source  is  neither  on  the  axis  nor  on  the  caustic  it  follows  that 
If  we  restrict  ourselves  to  the  axis  we  get 


^  •  \  "fTZ   I 


°°       ik(R*r) 

dV^    . 


>/Rr 


^a 


If  we  choose  ?      >    — -    we  have  a  point  of  stationary  phase  so  that 

O  fl 


,ik(RVr)  y?  [.K  (^  .  ^  )  (^/i4;  -  ^  ,")]' 
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The  behavior  of  u  with  respect  to  k  on  this  part  of  the  axis  is  in  no  way- 
distinguished  from  the  rest  of  the  space.  But  at  the  cusp  of  the  caustic  u  does 

have  a  singularity. 

2 
For  €q  <  >7  /*  there  is  no  stationary  point  and,  as  before,  we  replace  the 

variable  of  integration  ^  by  x/k  and  discard  all  terms  which  vanish  rapidly  with 

respect  to  k.  Then 

.ik(R+F)   /    ^^a 

"»^(>[a-  ^o^ 


u  -  2i  e-"       ,  ^ 


Thus  the  behavior  of  u  on  the  axis  above  and  below  the  cusp  of  the  caustic  is 
essentially  the   same,  nor  does  it  behave  very  differently  in  the  rest  of  tte 
space,  except  that  at  the  cusp  a  singularity  similar  to  a  line   source,  in  an 
asymptotic  sense,  seems  to  exist. 

On  the  axis.  In  the  iirjnediate  neighborhood  of  the  cusp,  the  function  R  +  r 
can  be  expanded  into 


where 


R+r     «     a+P€+YC 


a    -     a  +   ^^/a     ♦  2  v^^, 
^  a  -  a  gp 


Y     -     2v^^ 

2 

and  ^  Z  ^    /a.     Then 
o      V  a 


«o 


-(V^      (^iT^-l 


u     -     

8  Aexp["ik(a-  |^)J 


exp[ik(a+pC+YC^)]  d  /?  , 


.(v,..a,(M-^^     ^C,,(6)-13,,(M>. 
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where  C..  (6)   and  S..  (6)  represent  Hardy's  generadized  Airy  integral'-  J .     At  the 


cusp  the  above  reduces  to 


r(^)e^"/8 


8a 


-.i/U 


_k(i7  +ar_ 


r     ^2-, 
exp  kk(  yS"  ♦  -i  )  J 


-lA 


which  remainp  finite,  but  is  on  the  order  of  k  '  with  respect  to  k. 

A  different  behavior  can  be  detected  when  the  cusp  is  approached  along  the 
caustic.  We  will  assume  in  the  following  that  the  point  with  the  coordinates 
(4  f^  )  lies  on  the  caustic.  Then 


u 


n^ 


a 


,00 


<-00 


ik(R+r) 
2 dV^ 


where 


R  -  /(?-  ^a  -  ^J'^  *  ^«  ^a  » 


r  -  /(?-  ^a"  ^o  M  ^^  *  ''^'la*  ^^o"?  "  ^^^^o^^^^ 


The  function  R  +  r  has  two  points  at  which  it  is  stationary,  but  at  only  one  of 
these  do  both  the  first  and  second  derivatives  vanish.  Therefore  only  this 
point  will  be  considered.  Then  near  that  point  we  can  write 


R 


+  r  -  R"+  F  ♦  ^-^  I  ^  -  A^  ♦  ...  , 


where 


43  .  9^(R>r) 

A   "   —  '*■■  ■ 
34^ 


I. 


z 


BO  that 
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u 


2 
n 

'    R  r  4 

loo 

exp  ik 

6     J 

d€ 

2| 

r(i)   /^ 

"      /3Ri 

"4 

6^3 

• 

Near  the  cusp  of  the  caustic  this  expression  reduces  to 

2r(|)  r      .  11/6 


u 


It  may  be  observed  that  for  fixed  k,  u  Increases  as  the  cusp  Is  approached 
along  the  axis.  On  the  caustic  u  behaves  like  C*   ,  but  on  the  axis  u  behaves 
like  I""   •  However  at  a  fixed  point  not  on  the  caustic  u  is  of  the  order  k.         , 
but  of  the  order  k"  ' ^  on  the  caustic  away  from  the  cusp,  and  of  the  order  k"  ' 
at  the  cusp. 

In  conclusion  It  should  be  pointed  out  that  in  the  preceding  cases  only 
examples  of  boundary- value  problems  were  treated  for  which  the  Incoming  and  out- 
going waves  cancelled  on  the  surface.  The  cases  where  the  normal  derivatives 
cancel  could  be  treated  in  a  similar  manner  since  the  asymptotic  forms  of  the 
Green's  functions  G^  have  been  derived  for  both  the  paraboloid  of  revolution 
as  well  as  the  parabolic  cylinder.  The  resulting  integrations  reduce  to  the 
sane  evaluations  by  stationary  phase  techniques,  so  that  in  many  instances  the 
answers  could  almost  be  written  at  sight  upon  comparison  with  the  equivalent 
boundary-value  problem  of  the  first  kind. 
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